Internal mode mechanism for collective energy transport in extended systems 
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We study directed energy transport in homogeneous nonlinear extended systems in the presence of 
homogeneous ac forces and dissipation. We show that the mechanism responsible for unidirectional 
motion of topological excitations is the coupling of their internal and translation degrees of freedom. 
Our results lead to a selection rule for the existence of such motion based on resonances that explains 
earlier symmetry analysis of this phenomenon. The direction of motion is found to depend both on 
the initial and the relative phases of the two harmonic drivings, even in the presence of noise. 
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One intriguing phenomenon that is receiving much at- 
tention recently is net directed motion induced by zero 
average forces. Originally motivated by stochastic mod- 
els of biomolecular fbrownian) motors , deterministic 
ratchetlike systems 01 are being intensively studied, 
chiefly because of their many potential technological ap- 
plications 0. Many such models consist of one or two 
particles on a periodic, asymmetric potential and a peri- 
odic force (rocking ratchet Q). Later, the investigation 
was generalized to systems with many interacting parti- 
cles, from noisy soliton-bearing systems 0, to other 
spatially extended (stochastic and deterministic, over- 
damped and underdamped) systems, both theoretically 
and from a more applied 8] viewpoint. 

Among this class of problems, net transport in homoge- 
neous extended systems driven by homogeneous ac forces 
is particularly interesting. A paradigmatic example is the 
ac driven, damped sine-Gordon (sG) equation: 



4>tt - 4>xx + sin(0) = -f3(l>t + /(*)■ 



(1) 



A symmetry analysis, proposed for one-particle systems 
in (3 and extended to this problem 0, UiJ, indicated 
that a directed energy current appeared if f(t) broke 
the symmetry f(t) = -f(t + T/2), T being the pe- 
riod of the external driving. One such choice is f(t) = 
ei sm(St+S )+e 2 sm(mSt+5 +6) (H03 with S = tt/2), 
a case for which numerical simulations of the sG equa- 
tion confirmed the symmetry analysis results. In what 
follows, we will refer to So as the initial phase and to 
as the relative phase. Transport required a nonzero 
topological charge, implying the existence of sG solitons 
(kinks) in the system. In this respect, we stress that kink- 
mediated transport is impossible with only one harmonic 
for any value of the damping coefficient (3 [ljj . It was ar- 
gued in 01 that the observed rectification arises from 
the nonadiabatic excitation of internal kink modes and 
their interaction with the translational kink motion. This 
conjecture had no rigorous support; rather, it was based 



on plots of sG soliton evolution and on the failure of a 
collective coordinate (CC) approach with one degree 
of freedom, which assumed that sG solitons behave like 
rigid particles. An attempt to include the width degree 
of freedom has been recently presented in [lflj . where it 
was concluded that the directed energy current vanishes 
unless the width of the kink, l(t), is a dynamical variable. 
However, this condition is only a necessary one: l(t) is a 
dynamical variable in the one-harmonic case but the kink 
velocity is zero for any value of the damping as already 
mentioned Another point not accounted for in 
is the connection between the internal mode mechanism 
and the symmetry analysis, which also prohibits motion 
in other cases where l(t) is a dynamical variable. There- 
fore, the reasons for the phenomena observed in [jj 0] 
remained largely obscure. 

In this letter, a different CC approach allows us to 
identify the mechanism through which the width oscilla- 
tion drives the kink and its relation with the symmetry 
conditions. Furthermore, our theory predicts, and nu- 
merical simulations of Eq. (JJJ confirm, that the direction 
of motion depends on the initial phase of the driving, 
even in the presence of additive noise. Our CC theory is 
based on an Ansatz, proposed in |l4|. for the perturbed 
kink depending on two CC, X(t) and lit) (respectively, 
position and width of the kink). It is not difficult to show 
LLdl I If* Li! that the dynamics of these two CC is given 

by 



dP 
~dt 



-PP-qf(t), 



I 2 - 211 - 2(311 = Q 2 R l 2 



Mi 



a 



(2) 
(3) 



where the momentum P(t) = MqIqX /l(t), fin — 
l/(^/alo) with a = ir 2 /12 is the so-called Rice's fre- 
quency, and Mo = 8, q = 2n and Iq = 1 are, respectively, 
the dimensionless kink mass, topological charge and un- 
perturbed width. Equation (J2J can be solved exactly, 
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and in the large time limit (t 3> (3 1 ) yields 

P(t) — —y/e[ai sm(5t+So- X i)+a 2 sm(m5t+S + 9— X2 )], 

where e is merely a rescaling parameter in the 
perturbation expansion, to be determined later; 
Xi = arctan (5/(3) , X2 = arctan (m6/(3) , ai = 
qti/ yje((3 2 + S 2 ), and a 2 = qe 2 / ^Je((3 2 + m 2 5 2 ). As we 
are interested in the damped ((3 ^ 0) case and Eq. © 



cannot be solved in that case 0, , we will study it by 
a perturbative expansion, l(t) = Z + eZi(t) + ( 2 l 2 (t) +•■• ■ 
At order O(e), we obtain 

i'i(t) + 0h(t) + n 2 R h(t) = -n 2 R P 2 {t)l /2eMl (4) 

The key point is that, by substituting the expression of 
P(t) into Q, we see that the equation for li(t) contains 
harmonics of frequencies 2S, 2md and (m ± 1)6, i.e., 



h(t) + (3h(t) + fl 2 R h(t) = Ai+A 2 cos(26t + 26 - 2 X i) + A 3 cos(2mSt + 25 + 29 - 2 X2 ) 

+ Ai cos[(m - l)5t + 6-(x2- Xi)] - A ± cos[(m + l)6t + 26 + 9 - ( X2 +Xi)], 

where A\ = —A 2 — A 3 , and A 2 = Q,Ra\/A\fotMQ, A 3 = VlRol/A^faM^, A4 = —Cl Raids /2y/aM§. After transients 
elapse, we find 



hit) 



A 1 | A 2 sin(2(5t + 2^ -2xi + ^2) | A 3 sm{2mSt + 25p + 26 - 2 X 2 + ggm) 
V(^- 4(5 2 ) 2 + 4/3 2 (5 2 v/(fi| - 4m 2 (5 2 ) 2 + 4m 2 /? 2 <5 2 

A 4 sin[(m - 1)5* + - ( X2 - Xl ) + 6> m _i] A 4 sin[(m + l)£f + 2S Q + 9 - ( X2 + Xl ) + 9 m+1 ] 



n 2 



(5) 



y/(n% - (to - 1) 2 6 2 ) 2 + /3 2 (to - 1) 2 <5 2 



^(ft 2 , - (m + 1) 2 <5 2 ) 2 + /3 2 (m + 1) 2 S 2 



where 9 m = arctan — m 2 S 2 ) /mf3S\ . A cumbersome 
but otherwise trivial calculation yields the harmonics 
contained in l 2 (t), collected in Table [I] 

Next, we need to compute the average velocity over 
one period T — 2ir/5: In the CC approach, we use the 
definition of the momentum and find 



(*(*)> = i 



M l 



dt. 



(6) 



At O(e ), the averages (P(t)) and (Xo(t)) vanish triv- 
ially; therefore, net kink motion can only arise in next 
order. By straightforward calculations from Eqs. J5J and 
© we find for m = 2 that, for large enough times, 



„3r>2 



8M 3 (/? 2 + S 2 )^J(3 2 + 4S 2 



2cos[S Q -d+( X2 -2 Xl )-9 1 
^/{rf R -S 2 ) 2 +/3 2 S 2 



r 




(7) 



From Eq. (7) we see that for e to be small the prefactor 
on the rhs has to be much smaller than 1. A definite, 
verifiable prediction from this asymptotic expression is 
the existence of a nonzero velocity for m = 2, with a 
sinusoidal dependence on So and 9. This means that 
the velocity depends on both the initial and the relative 
phases; indeed, by letting So = Sto in Eq. (JIJ and chang- 
ing variables to t' = t + to, it can be immediately seen 
that an initial phase Sq is equivalent to a relative phase 
9' = 9 — (to — l)6o for a kink with its center shifted to 
xo + Vto- The dependence of the velocity on 9 agrees 
with (and explains) 0, ll(| , whereas the dependence on 



So is a totally new result. Nevertheless, these analytical 
results as well as the numerical simulations we present 
below strongly support the present conclusion. 

For the case m = 3, the average velocity is zero at all 
orders, a result confirmed by direct numerical simulation 
of the full sG equation (JIJ as we will see below. The rea- 
son can be understood by looking at Tabled For m = 3, 
the frequencies of the ac force (or the momentum) are 
odd harmonics (S and 36) , whereas the width of the kink 
oscillates only with even harmonics (2n5, n € N). This 
leads us to our main conclusion, namely the mechanism 
for the appearance of net motion and the corresponding 
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TABLE I: Harmonic content of the first contributions to the 
perturbative expansion of l(t). 



harmonic 


h 


h 


m 


28, 2m8, (m±l)5 


28, 48, 2mS, 4mS, (m ± 1)8, 
2(m ± 1)5, (m ± 3)8, (3m ± 1)8 


2 


8, 28, 38, 48 


8, 28, 38, 48, 58, 68, 78, 88 


3 


28, 48, 68 


28, 48, 68, 88, 108, 128 
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FIG. 1: Dependence of the kink velocity on the initial phase. 
Parameters are ei — £2 = 0.2, /3 = 0.05, 8 = 0.1. Relative 
phase 8 = n/2: solid line, CC theory; filled circles, simula- 
tion results. Relative phase 9 — 0: dashed line, CC theory; 
squares, simulation results. 



selection rules. Equations @ and © show that the force 
acts on the kink width through P 2 {t), whereas P{t) itself 
is in turn inversely proportional to l(t). This coupling is 
the responsible for the net kink motion, but for it to be 
actually possible, the harmonic content of the effective 
force P 2 {t) acting on the width degree of freedom must be 
able to resonate with it. This is evident from Eq. JfjJ, in 
which the integral is nonzero only if l(t) contains at least 
one of the harmonics of P(t). It is important to realize 
that this condition is much more restrictive than that 
found in [lflj . where only the necessity of l(t) being a dy- 
namic variable was pointed out. We have just seen that 
this is indeed necessary, but that additional, crucial res- 
onance conditions have to be fulfilled. Interestingly, our 
theory shows also that dissipation can change or even re- 
vert the kink velocity, see Eq. Q , in agreement with the 
numerical results in |fj, |l(J . A more detailed discussion 
of this point is forthcoming ■ 

These predictions from the CC approximation must 
be confirmed by numerical solution of the full partial 
differential equation This we do by using the 

Strauss- Vazquez scheme on systems of length of 

L = 100, 1000, with steps At = 0.01, Ax = 0.1, free 
boundary conditions, and a kink at rest as initial con- 
dition. Instead of the perturbative expressions (which 
are only qualitatively correct unless e << 1), to assess 
the validity of our theory we numerically integrate Eq. 
@ and the equation for the velocity obtained from the 
expression of P(t), Eq. (J2J), with a fourth-order Runge- 
Kutta method. Our main results are shown in Figs. ^421 
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FIG. 2: Dependence of the kink velocity on the initial phase 
for relative phase 9 — tt/2 in the deterministic (D — 0, empty 
circles) and the stochastic (D = 0.03, diamonds) cases. Other 
parameters are as in Fig. 



they fully confirm the accuracy, even quantitative, of our 
approach. Figure exhibits clearly the sinusoidal depen- 
dence of the velocity as a function of the initial phase. 
The dependence on 9 is also seen as a simple shift when 
changing from 9 = to 9 = tt/2. The agreement with the 
CC results is perfect. As a further check of the robustness 
of this dependence, following |i| we have simulated Eq. 
(|T|) with an additional additive gaussian white noise term 
with variance D. While one could in principle think that 
this noise would suppress the initial phase dependence, 
Fig. shows that the opposite is the case: the noise en- 
hances the dependence on the initial phase, increasing the 
maximum values of the velocity while keeping the same 
general sinusoidal dependence and the location of the ze- 
ros. It is tempting to conclude from this plot that the 
noise, at least if it is not very large (D « 1), assists 
the process of energy transfer between the width and the 
translation degrees of freedom, activating it. Finally, Fig. 
13 makes it clear that our main result, namely the inter- 
pretation of the physics of the problem, is indeed true, by 
showing the harmonic content of l(t) for m = 2 and 3. In 
this case, the agreement between our CC theory and the 
full numerical simulation of Eq. is indeed impressive, 
and validates firmly our resonance criterion for net kink 
motion. It is important to stress that the present th eory 
does not apply to the net motion found for m — 3 in [lfj . 
We have confirmed their result in our simulations, which 
allowed us to realize that this is an altogether different 
phenomenon: first, it appears only above a (moderately 
large, > 0.4) threshold amplitude, and second, it is 
induced by the kink wings, which are highly distorted in 
the process yielding the CC picture inappropriate (even 
kink-antikink pairs are created). 

In conclusion, we have found that the symmetry con- 
ditions set forth in 0, have their physical origin in 
the mechanism of the directed motion: the indirect ac- 
tion of the force through the coupling of the translational 
and width degrees of freedom. To make net motion pos- 
sible, this indirect driving has to resonate with the avail- 
able frequencies for the width. This interpretation does 
not contradict the non-existence of internal modes in sG 
kinks, shown in [l6|. because external forces can induce, 
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FIG. 3: Discrete Fourier Transform of the kink width. Upper 
panel: m = 2; lower panel: m = 3. Solid line: amplitude 
measured in simulations. Dashed line: numerical integration 
of the CC equations. Parameters are as in Fig. for relative 
phase 6 = tt/2 and So = —2.5. 



via excitation of certain phonons, behavior similar to the 
one expected from an intrinsic internal mode |l9l |20| . 
The fact that a force with only one harmonic would not 
drive the damped sG kink and that two harmonics 
are needed to simultaneously excite the width oscilla- 
tions and induce net motion fits nicely in this picture. 
On the other hand, this point rises the question as to the 
generality of our results, in view of the fact that most 
kink-bearing systems do have internal modes. To answer 
this question, we have studied the same problem in the 
framework of the </> 4 model, reaching the same conclu- 
sions Indeed, the intrinsic internal mode of 4 kinks 
makes the phenomenon even more noticeable, making us 
confident on the wide applicability of this work. 

Another important conclusion is the dependence of the 
velocity on the initial phase 8q, not mentioned in earlier 
work 0,^3- We note that this dependence allows much 
more flexibility in controlling the kink velocity, providing 
an alternative to the use of the relative phase suggested 
earlier. On the other hand, this may have important con- 
sequences for applications as a way of separating, e.g., 
fluxons in long Josephson junctions ||- Interestingly, 
such superconducting devices provide the best possible 
laboratory to verify our results. This experimental con- 
firmation is crucial in order to ascertain their applica- 
bility. Given the accuracy with which the sG equation 
describes long Josephson junctions, and the fact that an 
external force like the one proposed in this and earlier 
works [EE! is easy to implement, we hope that the cor- 
responding measurements will soon be carried out. A 
conclusive, positive verification of our theory would yield 
the picture we provide here very useful in that and re- 
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